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ABSTRACT 


It is shown that the influences of the thermal and tidal effects on Mercury's 
libration are in equilibrium with the periods of rotation and revolution of 
Mercury locked in the 3:2 resonant state. The suggestion by Liu that the solar 
gravitational couple on the thermal bulges accelerates Mercury's rotation is 
investigated and the production of mechanical energy to balance the dissipation 
of the bodily tides is discussed. It is possible for Mercury to rotate with two 
bulges as a solar thermal engine; the tidal effect causes this engine to function 
and its maximum power is close to 10 16 ergs per second. 
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THERMAL AND TIDAL EFFECT ON 


THE ROTATION OF MERCURY 


I. INTRODUCTION 

Liu (Liu, 1969; Liu, 1968) has suggested that the trapping of Mercury's 
rotational period into a 3:2 resonance lock with its orbital period was originally 
affected by the thermal contraction of the figure of Mercury during solidification. 
In the present paper the analysis of the contribution of the two thermal bulges to 
the dynamic stabilization of the planet's rotation is given. Attention is focused 
on the balance of the influence of the tidal and thermal effect on the libration of 
Mercury after solidification. 


II. BASIC EQUATIONS 

After solidification, the thermal contraction of the figure of Mercury from 
loss of heat must be exceedingly small and the thermal bulges considered by Liu 
(Liu, 1969; Liu, 1968) can grow because the apparent circulational motion at suc- 
cessive perihelia has been converted to a librational motion. The thermoelastic 
variation with time of the fractional difference in Mercury's equatorial moments 
of inertia about the state of hydrostatic equilibrium is 


B(t) - A(t) 
C(t) 


B(h) - A(h) 3(2 ~ 5/3 3 + 3/3 5 ) 
C ( h ) 100(1 + /3 + /3 2 ) 


• a - AT 


( 1 ) 


in which a is the coefficient of linear thermal expansion, AT the difference in 
surface temperature between the regions around the perihelion and the aphelion 
axes of Mercury, and /3 = (fRm - y)/$m where 5lm is Mercury's radius and y 
the depth of solar heating. The value of [B (h) - A (h)] /C (h) has been estimated 
as 3 x 10" 6 and the rate of increasing of [B(t) - A(t)]/C(t) is about 10" 21 
Sec. -1 (Liu and O'Keefe, 1965; Liu, 1968; Liu, 1969). Since the local tempera- 
ture regimes depend strongly on the orientation and the orbital position of 
Mercury relative to the Stun, a temperature phase effect is apparent with a phase 
lead of 35.5° (Liu, 1966b). Consequently, the angle between the axes of A(h) and 
the two thermal bulges is about 0.6°. 

In considering the bodily tidal torque we have recourse to the estimation by 
Jeffreys (Jeffreys, 1959). 


1 



N 


Km 6 ( 1 + e cos f ) 6 
a 6 ( 1 - e) 


• Sin 2 e 


(2) 


18 

5 



S 

k7rGp — • 

M 2 


• C 


(h) 


where G is the gravitational constant, p density, Ms mass of the Sun, Mm mass 
of Mercury, f the true anomaly, a the semimajor axis, e the orbital eccen- 
tricity, e the phase lag of the conventional equilibrium tides and k the Love 
number. 

The orientation of Mercury relative to the Sun, cp , is then governed by (Liu 
and O'Keefe, 1965) 


d 2 cp 

df 2 


2e sin f / dcp \ + 

1 + e cos f^df + ^J + 2(l+e cos f ) ^ 



( 1 - e 2 ) 3 
( 1 + e cos f ) 4 


(3) 


where 


_ B (h) - A(h) 
( h > C(h) 


AX 



• a - AT 


and n is the mean orbital angular velocity. 


III. ANALYSIS 

Let us choose the time r reckoned from perihelion and related to the 
orbital period divided by 2n 


T 


(O 


2 tan " 1 


/ 1 ~ e f e (1 - e 2 )^ Sin f 

\ 1 + e / an 2 1 + e cos f 


(4) 


As a new unknown function, we take the angle of rotation y between the axis 
of A(h) and the radius vector of the perihelion 


2 



f + cp 


(5) 


¥ - 


Equation (3) takes the form 

d 2 vj/ 3 (\h) + A M 
dr 2 2 


( 1 + e cos f ) 3 
(1 - e 2 ) 3 


Sin 2(\|/ - 


O 


N 


n 2 C 


(h) 


( 6 ) 


The solution of Equation (6) may be sought in the form 

— CIt + 8 


where ft is a constant and § is an unknown function. The resonances occur at 
2ft = m if m is an integer and S is, then, the angle of libration. To obtain an 
approximate expression for the libration of Mercury we may average it over a 
period of 2rr for m = 3. If the imperfection of elasticity adopted for the Earth 
(Jeffreys, 1959) is supposed to hold right up to the surface of Mercury, Equation 
(6) becomes 


d 2 S 2 • 1 A . ... 

+ ^ — (A.., . + AX) Sin 28 - - 10 11 

dr 2 2 (h) ' 


In the derivation of Equation (7), the following values were adopted: 


(7) 


e = 

0.2 





k = 

0.02 




2e = 

0.005 




G = 

6.7 

X 

o 

i 

00 

dyn. 

Cm 2 • g 2 

P = 

5.0 

g 

• Cm' 

-3 


n = 

1.2 

X 

10" 6 

rad. 

• Sec. -1 

M s/Mm = 

6.0 

X 

10 6 



fRm/a = 

4.1 

X 

10" 5 
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The average process also included the reversal of the tidal torque in direction 
during perihelion passage (Peale and Gold, 1965; Liu, 1966b). 

The first integral of Equation (7) is 


1 

2 



2 • 1 
4 


(X.^ h j + AX) cos 28 


10" 11 • 8 + E 0 


(8) 


where E 0 is a constant. 

Differentiating Equation (8) with respect to t , the result is 


_d_ _l_/d8\ 
dt 2 \dt/ 


2 - 1 


X 


(h) 


Gbs 28 


(9) 


2 - 1 


-n 


d(AX) 

dt 


., , d8 

Cos 28 - 10" u 1 n 2 • - 


2 • 1 


n 2 • AX • Sin 28 


dS 

dt 


The first and second term in the left side of Equation (9) correspond to kinetic 
and potential energies, respectively, of the system. In order to maintain the 
prescribed librational motion, energy must be added to the system because of 
tidal dissipation. Earlier investigations (Goldreich and Peale, 1968; Goldreich 
and Peale, 1966a; Goldreich and Peale, 1966b; Colombo, 1965; Colombo and 
Shapiro, 1966; Jefferys, 1966; Laslett and Sessler, 1966; Liu, 1966a; Liu and 
O'Keefe, 1965; Lutze and Abbitt, 1969) have shown that a sufficient deviation 
from axial symmetry would stabilize the libration of Mercury at this resonance 
state. Therefore the instantaneous deformation of the figure of Mercury into a 
hydrostatic equilibrium shape with X(h) = [B(h) - A(h)]/C(h) = 3 x 10“ 6 (Liu, 
1969) is sufficient to stabilize this resonance against the disruptive influence of 
the tidal or any other effects. However, the dynamic behavior and the evolution 
of the tidal or thermal process near this resonance condition are still unknown. 
To explain the physical feature of the locked-in motion, Liu (Liu, 1969; Liu, 

1968) has proposed that two bulges which arise on Mercury's surface from solar 
heating after the rotation becomes commensurable would increase the strength 
of the resonance and tend to prevent drifting out. This fact can be seen by com- 
puting the rate of working of the thermal and tidal effects. Equation (9) can be 
written in the form, 


4 



dE /dE 


dE\ 


dt 'thermal expansion ^ 'tidal effect * ^thermal bul 


dE 


(10) 


ges 


where 


dE 

dt 


thermal expansion 


2- 1 2 r _ d(A\) 

4 n ' 0)525 '“dr' c (t) 


/ dE\ 

'^^'tidal effect 


- 10" 11 • n 2 • — • C 
iU dt (O 


ME\ 

(dt) 


thermal bulges 


2 • 1 dS 

— 2 — ’ n 2 • AX. • Sin 2 S • — j-£- • t ^ 


For small angle of libration S, d(Ak)/dt = 10 21 Sec. -1 , and d8/dt < 10 -3 n, 
we obtain 


dE \ 

dt/ 


thermal expansion 


0(10 10 ergs. Sec. -1 ) 


dE 

dt 


t idal effect 


0(10 10 ergs. Sec. -1 ) 


The thermal expansion effect works against the tidal effect on Mercury's 
libration at a rate of 10 10 ergs. Sec -1 . It is of the same order of magnitude as 
the rate of the librational dissipation for the bodily tidal friction of Mercury. 

This result does not exclude entirely the possibility of equilibrium between the 
thermal expansion and the tidal effects but it must be admitted that the impres- 
sion remains that the tidal effect is the larger and that the rotation of Mercury 
is being retarded. For this reason, one is still inclined to be dissatisfied with 
the fact that the Liu-O’Keefe theory depends on a chance. This misunderstanding, 
however, can be easily cleared up. The most important feature of the rotational 
lock for the present case is that the main component is the one corresponding to 
the mechanical work done by the couple of thermal bulges. The accelerating 
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couple exerted by solar attraction of Mercury’s thermal bulges has been men- 
tioned. As Mercury seems to be the same age as the solar system, about 4.5 
billion years old, these bulges may well have existed for a large fraction of this 
time. Thus the value of AX at present time is about 10" 5 . If the deviation of 
the axis of the two thermal bulges from line of centers at positions f = 0° and 
± 45° increased 0.6° in the clockwise direction as seen from the north, the rate 
of working by the thermal bulges is 


(-) 

thermal 


bulges 


0(10 14 ergs. Sec." 1 ) 


which is about four orders of magnitude larger than that of the bodily tidal effect. 
As the amplitude of this deviation increases due to tidal effect, it will exert a 
couple which becomes progressively more effective in reducing the retardation 
until the thermal and tidal effects are in equilibrium with the orbital and rota- 
tional periods locked at the 3:2 resonant state. For an offset of the thermal 
bulges of 7r/4 from the radius vector at perihelion, this rate of working should 
multiply a factor of 100. Therefore Mercury may rotate as a celestial thermal 
engine; the tidal effect causes the engine to function and its maximum power is 
close to 10 16 ergs. Sec" 1 . It should be noted that the efficiency required of the 
planet Mercury as a thermal engine is very small. The maximum rate of working 
by the thermal bulges is only the small fraction 10" 7 of the total rate at which ef- 
fective heat energy is received from the Sun. The efficiency of an ideal thermal 
engine working on a temperature drop of 1 deg. C at 300 deg. K is about 10" 3 so 
there is a considerable margin to allow for irreversible processes or heat loss. 


IV. CONCLUDING REMARKS 

We have shown that the influences of the thermal and tidal effect on Mercury's 
libration are in equilibrium with the periods of rotation and revolution of Mercury 
locked in the 3:2 state. With regard to the interaction of the thermal effect, tidal 
friction and gravitation, it is concluded that the two thermal bulges on Mercury’s 
surface contribute dynamically to the stabilization of the planet's rotation. 
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